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In order to study consequences of the differences between the ionic-to-neutral and neutral- 
to-ionic transitions in the one-dimensional extended Peierls-Hubbard model with alternating 
potentials for the TTF-CA complex, we introduce a double pulse of oscillating electric field in 
the time-dependent Schrodinger equation and vary the interval between the two pulses as well 
as their strengths. When the dimerized ionic phase is photoexcited, the interference effect is 
clearly observed owing to the coherence of charge density and lattice displacements. Namely, 
the two pulses constructively interfere with each other if the interval is a multiple of the period 
of the optical lattice vibration, while they destructively interfere if the interval is a half-odd 
integer times the period, in the processes toward the neutral phase. The interference is strong 
especially when the pulse is strong and short because the coherence is also strong. Meanwhile, 
when the neutral phase is photoexcited, the interference effect is almost invisible or weakly 
observed when the pulse is weak. The photoinduced lattice oscillations are incoherent due to 
random phases. The strength of the interference caused by a double pulse is a key quantity 
to distinguish the two transitions and to evaluate the coherence of charge density and lattice 
displacements. 

KEYWORDS: photoinduced phase transition, neutral-ionic transition, TTF-CA, charge-transfer 
complex, time-dependent Schrodinger equation, dynamics 



1. Introduction 

In paper 1,1) it is suggested that we can induce the dynamics of charge density and that 
of lattice displacements simultaneously or separately by choosing the character of the applied 
electric field. In this sense, the electron-lattice coherence can be controlled. The charge-lattice- 
coupled dynamics is regarded as coherent if the ionicity and the staggered lattice displacements 
evolve simultaneously. In this case, the phases of the staggered lattice displacements are almost 
spatially uniform. In paper 11,^) we show the difficulty of spontaneously developing charge 
transfer to the ionic state possibly because the electron-electron interaction is responsible for 
the charge gap. In the present paper, we show that the degree of coherence is measured by 
using a double pulse and observing the resultant interference. Thus a double pulse experiment 
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would be a useful probe to characterize the photoinduced dynamics. 

In some cases, the degree of coherence may be guessed by comparing experimental data on 
different space and time scales, like photoreflectance with good space and time resolutions^^ 
and bulk-sensitive x-ray diffraction for the averaged structure'^) of the TTF-CA complex, al- 
though we have to pay attention to different experimental conditions including the energy of 
the pump light. For instance, in the spin-crossover complex [Fe(2pic)3]Cl2EtOH that shows 
a two-step thermal transition, although the intermediate-temperature phase is now known to 
have a symmetry-broken high-spin-low-spin configuration,^^ the x-ray diffraction does not so 
far show a corresponding photoinduced structural change as an averaged structure.®^ Such 
a difference between the thermo- and photo-induced structural changes may be related with 
different spatial correlations. In the TTF-CA complex, the correlation between the ionic- 
ity and the staggered lattice displacements or the spatial correlation of the latter would be 
important in dynamically controlling the transition in future. 

The coherence studied here is different from the coherence of the electronic excitations in 
quantum optics or from the quantum coherence needed for quantum computing. Even in pho- 
toinduced phase transitions, one from Mott insulator to metallic phases, for example, would 
proceed very fast without structural changes, so that it may have coherence of electronic 
origin. The ionic-to-ncutral transition is, however, necessarily accompanied with structural 
changes, so that there is basically no coherence of electronic origin. It might have some in- 
terference effect when the interval of a double pulse is comparable with the period of the 
excitonic oscillation, but it would be hardly related with the phase transition. As long as the 
phase transition is accompanied with structural changes, they determine the transition rate. 

Previously, we show a rapid and small-amplitude oscillation of the ionicity in the inset of 
Fig. 7 in ref. 8. However, only the slow and large-amplitude changes are accompanied with 
large structural changes. The rapid change corresponds to the excitonic oscillation, whose 
Fourier spectrum is rather broad [Fig. 8(a) of ref. 8]. Thus, the rapid oscillation cannot keep 
the electronic coherence on the time scale of the optical lattice vibration [Fig. 7(a) of ref. 8]. 
In short, the coherence studied here is basically of lattice origin, describing a state where 
the phases of the staggered lattice displacements are almost uniform in a wide spatial and 
temporal region inside the ionic domain. 

In paper I for the transition from the ionic phase, the ionicity and the staggered lattice 
displacements evolve on a common time scale when the pulse is strong and short, while they 
evolve on different time scales when the pulse is weak and long. If one looks at a certain 
site, when a strong and short pulse is applied, the lattice displacement oscillates without 
large friction before and after the neutral-ionic domain wall passes this site and the center of 
the oscillation is accordingly shifted. We use a double pulse in this paper to show that the 
interference effect survives for a wide range of pulse intervals when the split pulse is strong 
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and short, while it decays with increasing interval when the split pulse is weak and long. 
These facts have the same origin, i.e., the coherence is strong when the pulse is strong and 
short, while it is weak otherwise. In this sense, a double pulse experiment would be suitable 
for estimating the degree of coherence. In paper II for the transition from the neutral phase, it 
is shown that the polarization is hardly aligned in the ionic phase. We show here by a double 
pulse that the interference effect is much weaker than in the transition from the ionic phase. 
The coherence is then hardly achieved. In this way, the differences between the ionic-to-neutral 
transition in paper I and the neutral-to-ionic transition in paper II are characterized from the 
viewpoint of the coherence. 

In general, observing the coherence through the interference caused by a double pulse 
would be useful to study the relation between the dynamics of charge density and that of 
lattice displacements. If we find a difference between the simulated dynamics and the exper- 
imentally observed one through the interference effect, we would be able to construct better 
models. For instance, three-dimensionality and energy dissipation are missing in this series 
of papers, and their importance would be evaluated by comparing the numerical results and 
the experimental data. When we study in future how initially one-dimensional mctastable 
domains finally grow three-dimensionally, the interference effect may give a hint for clarifying 
the transition processes and for devising a scheme to control them. 

2. Extended Peierls-Hubbard Model with Alternating Potentials 

We use a one-dimensional extended Peierls-Hubbard model with alternating potentials at 
half filling, as in paper I, 

H = Hel + i^lat , (1) 



with 
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where, {ci,a) is the creation (annihilation) operator of a 7r-electron with spin a at site Z, 
'>^l,a = ^a^hf^i ''^l ~ '^',T '^'^i' i^ dimensionless lattice displacement of the Zth molecule 
along the chain from its equidistant position, and yi = ui+i — ui. The distance between the 
Zth and {I + l)th molecules is then given by n = ro(l -|- ui+i — ui), where ro is the averaged 
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distance between the neighboring molecules along the chain. The other notations are also the 
same as in paper I. We numerically solve the time-dependent Schrodinger equation with the 
help of the unrestricted Hartree-Fock approximation for the electronic part, and the classical 
equation of motion for the lattice part, as described in paper I. Note that random numbers 
are added to the initial yi and iii values according to the Boltzmann distribution at a fictitious 
temperature T. 

Photoexcitations are introduced by modifying the transfer integral with the Peierls phase. 
In contrast to papers I and II, a double pulse of electric field is used. The first pulse is given 

by 

^(0 = -^sinwext* , (4) 
for < t < NextText with integer Next, Text = 27r/a;ext, and the second pulse by 

^i't) = — 1^ sintJcxt(i - i2nd) , (5) 

for t2nd < t < t2nd + -^cxtT'cxt- E{t) IS zero otherwise. If t2nd is smaller than A^cxt^cxt; wc add 
up the two terms for t2nd < t < A'ext^cxt- Therefore, if t2nd is zero, it coincides with the single 
pulse E{t) = Eext sinwext* with amplitude Eext and frequency coext- 

3. Results and Discussions 

We use iV=100, io=0.17eV, C/=1.528eV, y=0.604eV (when we start from the ionic phase) 
or F=0.600eV (when we start from the neutral phase), d=2.716eV, /?2=8.54eV, A;i=4.86eV, 
/u2=3400eV, and the bare phonon energy Wopt = (l/ro)(2A;i/mr)^''^=0.0192eV, and impose 
the periodic boundary condition. Here is the reduced mass for the donor and accep- 
tor molecules. With these parameters, the dimerized ionic phase is stable (metastable) and 
the neutral phase is metastable (stable) for y=0.604eV (y=0.600eV). These parameters 
are the same as in ref. 8, so that the bare phonon energy Wopt used here is about five 
times higher than the optical phonon energy of the TTF-CA complex. The ionicity is de- 
fined as p = 1 + (1/A^) Z^ili(— !)'("/). The staggered lattice displacement is defined as 

yst = ii/N)Elii-iyyi- 

3.1 lonic-to-neutral transition 

In this subsection, we adopt the ionic state as the initial condition for the time-dependent 
Schrodinger equation. If we use sufficiently weak fields with a given pulse duration, we do not 
observe any substantial change in the ionicity, irrespective of whether the field is split into two 
or not, and no matter how long the interval is between the split pulses. If we use sufficiently 
strong fields on the other hand, the transition always takes place to the neutral phase. There 
is thus a range of the field strengths for which the evolution of the ionicity depends on the 
interval, as shown below. By adopting the field strength in this range, we vary the interval 
between the two pulses to show the time dependence of the ionicity in Figs. 1(a) and 1(c) and 
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Fig. 1. Time dependence of (a) the ionicity and (b) the staggered lattice displacement in case the 
transition is achieved, and that of (c) the ionicity and (d) the staggered lattice displacement 
in case the transition is not achieved. The interval a;opti2nd is varied between the two pulses of 
'^ext/wopt=28, eai^ext/<^opt=8, and A^ext=3. The initial state is ionic at T/t=10~^. 



that of the staggered lattice displacement in Figs. 1(b) and 1(d). As shown in Figs. 1(b) and 
1(d), the period of the optical lattice vibration Topt is slightly less than 4/a;opt in the ionic 
phase. 

When the interval between the two pulses t2nd is a multiple of Topt, the transition takes 
place [Figs. 1(a) and (b)]. Here, the field is so strong that the single pulse (t2nd=0) first 
reverses the polarization [Fig. 1(b)] and the reduced ionicity is partially restored after one 
period [Fig. 1(a)]. After two periods, the ionicity is again partially restored but its intensity 
is weaker than that after one period. The staggered lattice displacement shows decaying 
oscillation. For t2nd being a finite multiple of Topt, the overshoot of the staggered lattice 
displacement after the second pulse is quickly weakened and soon becomes invisible with 
increasing t2nd- Then, the staggered lattice displacement decays almost simultaneously with 
the ionicity. The speed of the transition slightly and gradually slows down with increasing 

i2nd- 
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On the other hand, when the interval is a half-odd integer times Topt, the transition 
is suppressed by the destructive interference [Figs. 1(c) and (d)]. After the second pulse is 
applied, the ionicity slightly drops, but it never reaches a value in the neutral phase [Fig. 1(c)]. 
In Fig. 1(d), one sees that the first and second pulses hit the staggered lattice displacement 
in the opposite directions. The first pulse decreases the staggered lattice displacement, while 
the second one increases it. This is why the interference is destructive. 

Here we used by chance, for the frequency of the electric field, a multiple of the bare 
phonon energy. The results are irrespective of whether their ratio is an integer or not. The 
actual frequency of the optical lattice vibration is different from (higher because of the k2 term 
than) the bare phonon energy, so that the frequency of the electric field is incommensurate 
with the frequency of the optical lattice vibration. 

When the ionic state is photocxcited, the final state is clearly identified as either the ionic 
state or the neutral state, i.e., no intermediate state is produced, as demonstrated in paper I. 
Here we show the ionicity at LJopti=50 in Fig. 2, as a function of the interval between the two 
pulses t2nd- The ionicity takes basically two values, as expected, about 0.8 if the final state is 
ionic and about 0.3 if neutral. The high or low value appears almost periodically in Fig. 2(a), 
and its period is about that of the optical lattice vibration. For i^opti2nd ^30 the ionicity in the 
neutral state becomes significantly larger than 0.3. This is simply because the transition is not 
completed yet. If we plot the ionicity, e.g., at u;opti=100 instead of a;opti=50, it is again about 
0.3. For short pulses such as we presently adopt (A'^oxt=3), the range of the field strengths 
for which the interference is clearly observed is wide. In Fig. 2(b), we use a stronger pulse 
to show the similar result. Although the final ionicity is not as periodic as in Fig. 2(a), the 
interference is clearly seen even for a;opti2nd —30. The deviation from the periodicity is due to 
the fact that the staggered lattice displacement is not a pure sinusoidal function with a single 
frequency (not shown). Because other frequencies are substantially involved for the stronger 
pulse, the amplitude of the oscillation of the staggered lattice displacement changes with time 
after the first pulse. When the amplitude is small, the interference effect is suppressed as in 
Fig. 2(b) for a;opti2nd —20. We have performed extensive numerical calculations with different 
strengths, durations and intervals of pulses. The periodicity as a function of the interval is lost 
as shown here in an accidental manner. The oscillation of the staggered lattice displacement 
contains different frequency components around the main one, which seem to have stochastic 
distribution depending on the random numbers added to the initial lattice variables. 

Now we use a weaker and longer pulse to show how the interference effect is weakened. The 
range of the field strengths for which the final state depends on the interval is substantially 
narrower than before. The time dependence of the ionicity is shown in Figs. 3(a) and 3(c), 
and that of the staggered lattice displacement in Figs. 3(b) and 3(d). After the first pulse, 
the oscillation of the staggered lattice displacement decays faster [Fig. 3(b)]. Even when the 
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Fig. 2. lonicity at a;optt=50, as a function of the interval between the two pulses tt;opt*2nd5 (a) for 
ea£'ext/'^opt=8, and (b) for ea£'ext/'^opt=10- The electric field with a;ext/'^opt=28 and A'ext=3 is 
applied to the ionic phase at T/t=10~^. The final state is either ionic or neutral. 



second pulse is applied only Topt after the first pulse, the staggered lattice displacement shows 
only a weak overshoot and its oscillation in the neutral phase is incoherent. When the second 
pulse is applied 2ropt after the first pulse, the transition is not smooth any more and it takes 
a much longer time for the staggered lattice displacement to disappear. The evolution of the 
ionicity after the second pulse is also rather irregular and much slower [Fig. 3(a)] than before. 
Nevertheless, the coherence is still clearly observed. When the interval is a half-odd integer 
times Topt, the second pulse little affects the ionicity [Fig. 3(c)] and the staggered lattice 
displacement [Fig. 3(d)]. 

With the same strength and duration of the pulse as in Fig. 3, we show the ionicity at 
a;opti=25 in Fig. 4(a), as a function of the interval between the two pulses t2nd- Again, the 
final ionicity takes basically two values, about 0.8 and 0.3. The rather large values around 0.5 
in the neutral state for a;opti2nd >10 is due to the fact that the transition is not completed 
yet. If we plot the really final ionicity instead, they become about 0.3. In any case, the final 
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Fig. 3. Time dependence of (a) the ionicity and (b) the staggered lattice displacement in case the 
transition is achieved, and that of (c) the ionicity and (d) the staggered lattice displacement 
in case the transition is not achieved. The interval Wopti2nd is varied between the two pulses of 

Wext/<j^opt=28, eai?cxt/wopt=6.2, and iVcxt=5, which are weaker and longer than in Fig. 1. The 
initial state is ionic at T/t=10~^. The coherence is slightly weaker than in Fig. 1. 
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ionicity largely deviates from the periodic function. This is again due to the fact that the 
staggered lattice displacement is not purely sinusoidal with a single frequency [Fig. 3(b)]. 
Because other frequencies are substantially involved for the longer pulse (even if it is weak), 
the amplitude of the oscillation of the staggered lattice displacement substantially changes 
with time after the first pulse, as clearly seen in Fig. 3(b). As the pulse becomes longer (and 
weaker to maintain the number of absorbed photons), the final ionicity more strongly deviates 
from the periodic function [Fig. 4(b)]. In addition, the range of the field strengths for which 
the final state depends on the interval becomes narrow. 

3.2 Neutral-to-ionic transition 

In this subsection, the initial state is neutral and its lattice displacements are distributed 
around zero in equilibrium. In contrast to the previous subsection, the clear interference effect 
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Fig. 4. lonicity at a;optt=25, as a function of the interval between the two pulses tt;opti2nd5 (a) for 
ea£'ext/'^opt=6-2 and A'ext=5, and (b) for eai^ext/'^opt=5-5 and Arext=10. The electric field with 
t^ext/t^opt=28 is applied to the ionic phase at T/t=lQ~^. The final state is still either ionic or 
neutral. 



is not seen with any field strength. The final lonicity is much more insensitive to the interval 
between the two pulses. Nonetheless, a small-amplitude oscillation is found as a function of the 
interval, as shown later. We show the time dependence of the ionicity in Fig. 5(a) and that of 
the staggered lattice displacement in Fig. 5(b) with different intervals. The first pulse is applied 
for < a^opti < -^ext'^optrext=2.1. Then, for a;opti2nd <2, the two pulses directly interfere with 
each other to show the superficial dependence on a;opti2nd- Otherwise, the ionicity is insensitive 
to a;opti2nd [Fig- 5(a)]. To assign a;opti2nd to each curve, refer to Fig. 6(a) later. The reason 
for the insensitivity is rather obvious when one sees the evolution of the staggered lattice 
displacement [Fig. 5(b)]. Since the initial lattice displacements are distributed around zero, 
they start to oscillate in a random manner. Some displacements contribute to positive yst, 
others to negative yst- Due to the cancellation, the magnitude of yst is small. The second pulse 
also affects the displacements in a very similar manner. Then, the interference effect is hardly 
seen when the staggered lattice displacements are spatially averaged. 
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Fig. 5. Time dependence of (a) the ionicity and (b) the staggered lattice displacement. The interval 
a;opti2nd is varied between the two pulses of a;ext/'^opt=30, ea£'ext/'^opt=3, and A^ext=10. The 
initial state is neutral at T/t=10~^. The coherence is almost completely lost. 



When the neutral state is photoexcited, the final state continuously changes with the 
strength and duration of the pulse, as demonstrated in paper II. The ionicity at LOoptt=100 is 
shown in Fig. 6, as a function of the interval between the two pulses t2nd- The final ionicity 
continuously changes with t2nd and also with the field strength, as expected. With the stronger 
field, it is larger on average [Fig. 6(b) compared with Fig. 6(a)]. With the weaker field, the 
final ionicity appears to oscillate [Fig. 6(a)], although the amplitude of the oscillation is 
much smaller than that in the previous subsection. An interference effect is experimentally 
observed in the photoinduced dynamics from the neutral phase,^) though it is not compared 
yet with the interference effect expected in the photoinduced dynamics from the ionic phase. 
The numerically observed period is about a half of Topt in contrast to Topt in the previous 
subsection. 

To demonstrate that this small-amplitude oscillation is generally found with weak fields, 
we use shorter pulses to show the ionicity at a;optt=100 in Fig. 7, as a function of the interval. 
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Fig. 6. lonicity at u)optt=WO, as a function of the interval between the two pulses Wopt*2nd) (a) for 
ea£'ext/'^opt=3, and (b) for ea£'ext/'^opt=5- The electric field with a;ext/'^opt=30 and A^ext=10 is 
applied to the neutral phase at T/t=10~^. The final state changes continuously with eai^ext/^^opt- 



The results are similar. Namely, the final ionicity continuously changes with f2nd a^nd with the 
field strength. The final ionicity is larger for the stronger field. Again with the weaker field, 
the final ionicity oscillates [Fig. 7(a)]. The amplitude and period of this oscillation is similar 
to those in Fig. 6(a). 

The origin of the small-amplitude oscillation of the spatially averaged, staggered lattice 
displacement is rather obvious when one locally sees the staggered lattice displacement. The 
space and time evolution of the ionicity and the staggered lattice displacement is shown 
in Fig. 8. The horizontal component of each bar represents the local ionicity pi defined as 
Pl = 1 + (—!)'(— (n/_i) + 2{ni) — (n/+i))/4. The vertical component gives the local staggered 
lattice displacement ygti defined as ygti = (— !)'(— y^-i + 2yi — y/+i)/4. The bars are shown 
on all sites in the first 40 sites. The oscillating electric field is applied for < Uoptt < 1.05 
and 11 < LOoptt < 12.05. It is clearly seen that the first pulse induces oscillations of the 
local staggered lattice displacements and that their phases are randomly distributed over 
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Fig. 7. lonicity at u)optt=WO, as a function of the interval between the two pulses Wopt*2nd) (a) for 
ea£'ext/'^opt=5, and (b) for eai^ext/'^opt=10- The electric field with a;ext/'^opt=30 and A'ext=5, 
which is stronger and shorter than in Fig. 6, is applied to the neutral phase at T/t=10~^. 



the system. Because nearby oscillations with different phases compete with each other, each 
oscillation is not regular if one sees it very carefully. Nevertheless, the local oscillations persist 
in this time scale. Therefore, if the second pulse is applied a multiple of Topt after the first 
pulse, as shown here, it enhances the local oscillations so that the (spatially averaged) final 
ionicity becomes a little bit larger [compare a;opti2nd=ll in Fig. 7(a) with the nearby points 
at a;opti2nd=10-5 and 11.5.]. 

It should be noted that, if the second pulse is applied a half-odd integer times Topt after the 
first pulse, the final ionicity becomes a little bit larger again owing to constructive interference, 
in contrast to the destructive interference when the ionic state is photoexcited. One can 
imagine that a pendulum starts to oscillate from the bottom of a potential by a pulse of 
external force. When the second pulse is applied a half-odd integer times Topt after the first 
pulse, the pendulum is at the bottom of the potential, but its velocity is opposite to the initial 
one. Note that the electric field induces charge transfer in both directions. The direction to 
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Fig. 8. Correlation between the staggered lattice displacement yst i (the vertical component of the 
bar) and the ionicity pj (the horizontal component of the bar) , as a function of the site index I 
and the elapsing time t multiplied by o^opt- The electric field with eai?oxt/^opt=5, A^oxt=5 and 
of frequency (A;ext/i^opt=30 is split into two with an interval of Wopt^2nd=ll and is applied to the 
neutral phase at T/t=10~^. 



which the pendulum is hit makes no sense. Thus, the two pulses constructively interfere with 
each other also when the interval is a half-odd integer times Topt- As a consequence, the 
small-amplitude oscillation appears in the (spatially averaged) final ionicity and its period is 
about a half of Topt [Figs. 6(a) and 7(a)]. If the electric field is substantially larger, however, 
the cancellation among different phases of local oscillations is more effective, so that the 
small-amplitude oscillation becomes invisible [Figs. 6(b) and 7(b)]. 

4. Conclusions 

In papers I and II, we have shown the qualitative differences between the ionic-to-neutral^^ 
and neutral-to-ionic^) transitions and the field-strength-dependent relative dynamics of charge 
density and lattice displacements^^ in the one-dimensional extended Peierls-Hubbard model 
with alternating potentials. To study their consequences and to help us devise a scheme 
to control the transition dynamics, we here use a double pulse of oscillating electric field 
with different intervals in solving the time-dependent Schrodinger equation for the mean-field 
electronic wave function. 

When the photoexcitation induces the transition from the dimerized ionic phase with 
ferroelectrically ordered polarizations to the neutral phase with almost equidistant molecules, 
the interference effect is clearly observed due to the coherence of the ionicity and the staggered 
lattice displacements. The transition is promoted when the interval between the two pulses 
is a multiple of the period of the optical lattice vibration, while it is suppressed when the 
interval is a half-odd integer times the period. The lattice vibrations around the dimerized 
structure induced by the first and second pulses are in phase or out of phase, respectively. 
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In paper I, the infrared light nearly resonating with the optical lattice vibrations is shown to 
induce the transition as well as the light for creating excitons. The coupling between charge 
density and lattice displacements is strong enough to transfer charge by lattice oscillation. In 
this paper, the relative phase between the lattice vibrations and the oscillating electric field 
is shown important for the photoinduccd transition. 

In paper I, the correlation between the dynamics of the charge density and that of the 
lattice displacements is shown to depend on the strength of the pulse. When the pulse is strong 
and short, the charge transfer takes place on the same time scale with the disappearance 
of dimerization. When the pulse is weak and long, the dimerization-induccd polarization is 
disordered much before the charge transfer. In the present calculations with a double pulse, 
the coherence of the ionicity and the staggered lattice displacements is also shown to depend 
on the strength of the pulse. When the pulse is strong and short, the coherence is so strong 
that the interference effect is observed even for the intervals much longer than the period of 
the optical lattice vibration. When the pulse is weak and long, the coherence is weakened, 
and then the interference effect fades away with increasing interval. These findings have the 
same origin. Therefore, the double-pulse experiment would be a good tool to investigate the 
coherence of the ionicity and the staggered lattice displacements. 

When the photoexcitation induces the transition from the neutral phase with thermally 
disordered lattice displacements to the ionic phase, a very weak interference effect is observed, 
although the staggered displacements are spatially incoherent. The induced ionic state is 
stabilized by staggered lattice displacements, but they need not be spatially coherent in one 
dimension. The induced polarizations are thus disordered. The final ionicity depends weakly 
on the interval between the two pulses. Nevertheless, when the field is not so strong, the 
weak interference effect is certainly observed. The final ionicity oscillates, as a function of 
the interval, with small amplitude and twice as frequently as observed in the ionic-to-neutral 
transition. To realize the constructive interference, the second pulse must be applied when 
the lattice displacements have the maximum velocity. It is hard to observe the interference if 
the staggered lattice displacements are spatially averaged, but it is relatively easy if they are 
measured locally. As the pulse is strengthened, different phases of nearby oscillations compete 
more strongly with each other and suppress the interference effect. 

When quantitatively comparing the numerical results with the experimental observation, 
three-dimensionality and energy dissipation are not negligible. The double pulse experiment 
may be very useful in that it allows comparisons from the viewpoint of coherence. Because 
the information on the coherence of the ionicity and the staggered lattice displacements is 
obtained by the double pulse experiment, the relevance of the isolated one-dimensional model 
system can be discussed in a straightforward manner. The relevance may depend on the time 
scale of observation. 



14/16 



J. Phys. Soc. Jpn. Full Paper 

Acknowledgement 

The author is grateful to H. Okamoto for showing his data prior to publication and 
for enhghtening discussions. This work was supported by Grants-in-Aid for Scientific Re- 
search (C) (No. 15540354), for Scientific Research on Priority Area "Molecular Conductors" 
(No. 15073224), for Creative Scientific Research (No. 15GS0216), and NAREGI Nanoscience 
Project from the Ministry of Education, Culture, Sports, Science and Technology, Japan. 



15/16 



J. Phys. Soc. Jpn. Full Paper 

References 

1) K. Yoncmitsu: J. Phys. Soc. Jpn. 73 (2004) former of the previous papers. 

2) K. Yonemitsu: J. Phys. Soc. Jpn. 73 (2004) latter of the previous papers. 

3) S. Iwai, S. Tanaka, K. Fujinuma, H. Kishida, H. Okamoto and Y. Tokura: Phys. Rev. Lett. 88 
(2002) 057402. 

4) L. Guerin, E. Collet, M.-H. Lemee-Cailleau, M. Buron-Le Cointe, H. Cailleau, A. Plech, M. Wulff, 
S. Koshihara and T. Luty: Chem. Phys. 299 (2004) 163. 

5) D. Chernyshov, M. Hostettler, K.-W. Tornroos and H.-B. Biirgi: Angew. Chem. Int. Ed. 42 (2003) 
3825. 

6) N. Huby, L. Guerin, E. Collet, L. Toupet, J.-C. Ameline, H. Cailleau, T. Roisnel, T. Tayagaki and 
K. Tanaka: Phys. Rev. B 69 (2004) 020101. 

7) T. Luty and K. Yonemitsu: J. Phys. Soc. Jpn. 73 (2004) 1237. 

8) N. Miyashita, M. Kuwabara and K. Yonemitsu: J. Phys. Soc. Jpn. 72 (2003) 2282. 

9) H. Okamoto: private communications. 



16/16 



